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Abstract 

We introduce an entanglement criterion to exclude full separability of quantum 
states. We present numerical evidence that the criterion is necessary and sufficient 
for the class of GHZ diagonal three-qubit states and estimate the volume of bound 
entangled states within this class. Finally, we extend our approach to bound 
entangled states which are not GHZ diagonal. 



1. Introduction 

Entanglement is believed to be an important resource in quantum information 
processing and consequently many works are devoted to its characterization [l|, 
|2j. This characterization becomes complicated, if more than two particles are 
entangled, since then different classes of entanglement exist. Some methods to 
distinguish between the different classes have been presented in the literature 0- 



17|. However, no general solution of the problem is known, not even for specific 
families of states. 

In this paper we present a criterion for the verification of multi-qubit entan- 
glement. The criterion is formulated as a set of simple inequalities for the matrix 
elements of the state. Our approach is inspired by some existing criteria 11211 . and 
also by some recent works on entanglement in the family of graph states [la. fl7|. 
We consider then density matrices which are diagonal in terms of Greenberger- 
Horne-Zeilinger (GHZ) states and investigate the optimality of our criterion. For 
that, we derive methods to prove that a given state is separable and it turns out 
that all GHZ diagonal states under scrutiny are either detected by our criterion 
or proven to be separable. This allows to estimate the volume of so-called bound 
entangled states in this class of states. Bound entanglement is a weak form of 
Preprint submitted to Elsevier January 4, 2011 



entanglement, where no pure state entanglement can be distilled from; and this 
phenomenon is central to many open problems in quantum information theory 
[I] . Finally, we discuss with the help of an example how our ideas can be used to 
characterize bound entanglement close to the three-qubit W state. 



2. Definitions and statement of the problem 

Before explaining our separability criteria, we introduce the notation and give 
some examples of existing separability criteria and interesting quantum states. 

2.1. Separability and entanglement 

We consider an iV-particle system with Hilbert space %tot = % . Any ma- 
trix g acting on %tot which is hermitian (g = g^), has no negative eigenvalues 
(g > 0) and is normalized [Tr(g) = 1] is a valid density matrix of some quantum 
state. By definition, a state is fully separable, if it can be written as a convex 
combination of product states, 

g = ^2pk\a k }(a k \ \b k }(b k \ ® ••• ® \z k )(z k \, (1) 
k 

where the p k are non-negative (p k > 0) and normalized (YlkPk = 1); m other 
words, they form a probability distribution. If a state cannot be written as in 
Eq. (pQ) it is entangled in some sense. In general, it is not easy to check whether 
a given quantum state is fully separable or not, see Refs. [3H13I] for some existing 



sufficient criteria for entanglement and Refs. [181 . |19( for numerical tests to prove 
separability. Especially if a state is only weakly entangled, proving entanglement 
is not straightforward. 

In this paper we will derive criteria, which allow to prove that a state is not 
fully separable and hence contains some entanglement. It should be stressed 
that there are more refined notions of entanglement for multiparticle systems 
(e.g., genuine multipartite entanglement) which can be more relevant for special 
situations (e.g., experiments). We will, however, concentrate on full separability 
as the basic definition, for a discussion of recent results on the other classifications 
see the reviews in Refs. Moreover, we will focus our discussion on the case 

of three qubits, but our results can directly be generalized to more particles. 

2.2. A criterion for full separability in terms of the matrix elements 

Let us explain a separability criterion for full separability derived in Ref . [l2j , 
which we will generalize later. For three qubits, consider the 8 x 8-density matrix 
with entries Qij. Here and in the following, we always order the basis vectors in the 
canonical way, {|000), |001), 1 010) , | 111)}. Then, if the state is fully separable, 
the entries fulfill 

1 01,8 1 < ^£2,2 £3,3 £4,4 £5,5 £6,6 £7.7- (2) 
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The idea of the proof is as follows: It is easy to check that for a pure fully 
separable state equality holds in Eq. t^. Then, \gifi\ is convex in g [a function 
f(g) is convex if f[pgi + (1 — p)g2\ < pf(gi) + (1 — p)f(g2)], while the right-hand 
side of Eq. ([2]) is concave [i.e., f[pgi + (1 - p)q2[ > pf(gi) + (1 - p)f(g2)]- This 
implies the bound for mixed fully separable states. 

The right-hand side of the inequality ([2]) may be replaced by other expres- 
sions, for instance, \q\ % %\ < -^^2,2^3,3^5,5^8,8 holds for separable states. Then, 
this criterion is able to detect the entanglement in a family of states which are 
separable with respect to any bipartition, but not fully separable [l^]. This is 
then a bound entangled state and the criterion improves existing criteria for this 
family significantly. 

In the following, we will extend the criterion in Eq. ([2]) by considering more 
than one off-diagonal element on the left-hand side. We will also discuss the 
optimal choice of the right-hand side. 



2.3. A bound entangled state 

Let us now introduce a bound entangled state which motivates our approach. 
Consider the family of three-qubit density matrices introduced by A. Kay 16|, 
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This matrix is a valid quantum state for a > 2 and it is separable for d > 
2V2. This has been proved in Ref. [161 ] by writing down an explicit separable 
decomposition as in Eq. (pQ). Furthermore, using the algorithm outlined in Ref. [H] 
it was shown numerically that the state is entangled for 2 < a < 2.828. Note that 
for any d > 2 the state has a positive partial transpose (PPT) for any bipartition. 
This does not only imply that no entanglement can be distilled from i10; moreover, 
since the state is diagonal in the GHZ basis (see also below), it also implies that 



At this point, one should mention that undistillable entanglement can arise in multiparticle 
systems in a simple way: If g is entangled with respect to one bipartition, but PPT with respect 
to another bipartition, this implies already that g is entangled, but multipartite undistillable. 
In this paper, however, we use a more restrictive definition: we consider bound entangled 
states, that are separable for any bipartition, but not fully separable. This means that they are 
undistillable, even if arbitrary parties join. 
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the state is separable for all bipartitions. Namely, from the fact that a two- 
qubit Bell diagonal state with a positive partial transpose is separable, one can 
conclude that a GHZ diagonal state that has a positive partial transposition is 
also separable for that partition [io| . 

The state in Eq. ([3]) is an example of a GHZ diagonal state. These are of the 
form 

8 

g = Y J Pk\GHZ k ){GHZ k \, (4) 

k=l 

where the GHZ state basis consists of the eight vectors \GHZ k ) = IOX2X3) ± 
I1X2X3) where Xj,Xj £ {0, 1} and Xj 7^ Xj. Alternatively, the GHZ diagonal states 
can be written as 

e = - [1 + \ 2 ZZ1 + X 3 Z1Z + X A tZZ + X 5 XXX + X&YYX + X 7 YXY + X 8 XYY] , 
8 

(5) 

where X, Y, and Z denote the the Pauli matrices and tensor product symbols have 
been omitted. The observables occurring in Eq. ([5]) are the so-called stabilizing 
observables of the GHZ states, see Ref. 



2l| for a further discussion. 



GHZ diag onal states have been intensively discussed in the literature before 
0, l3, 2(], 22] and are interesting from several perspectives: They have a simple 



structure, since only the diagonal and the anti-diagonal elements of the matrix 
can be nonzero, and they occur naturally in certain types of decoherence processes 
[23| . Furthermore, any state can be transformed to a GHZ diagonal state (without 
changing the fidelities of the GHZ states) by local operations 0, 21]. This means 



that if the remaining GHZ diagonal state is entangled, the initial state must have 
been entangled, too. On the other hand, if one has a GHZ diagonal state g± and 
a second separable state Q2 which is mapped to £1 by these local operations, then 
£1 must be separable, too. 

The state qak is not detected by the criterion in Eq. ([2]) or variants thereof. 
The reason lies in the fact that Eq. (|2|) considers only single off diagonal elements, 
and neglects the (phase) relations between them. It is the main purpose of 
this paper to develop an improvement of Eq. ([2]) which takes into account all 
off diagonal elements at the same time. This will finally prove analytically that 
the state qak is entangled if 2 < a < 2y/2. 

3. The separability criterion 

In this section we will now formulate the separability criterion. We will restrict 
our attention to three qubits, but it should be stressed that our approach can 
straightforwardly be generalized to an arbitrary number of qubits. 

To start, consider a pure product state, 

\<f>) = (ci|0) + ® (c 2 |0) + s 2 |l» ® (c 3 |0> + s 3 |l». (6) 
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Here, the complex coefficients q and Sj fulfill the normalization |cj| 2 + \si\ 2 = 1. 
Let us consider the corresponding 8x8 density matrix g = |0) (0| . The off diagonal 
elements can be written as 

2i,8 = c^s^s = ^ Q+6+c \ 

£2,7 = C1C2S3S1S2C3 = Ke M , 

£3,6 = C1S2C3S1C2S3 = Ke , 

£5,4 = siC2C 3 c 1 s 2 s 3 = Ke . (7) 

Here, we used the notations a = <j>\ + 02 — 03, b = <\>\ — 02 + 03, and c = 
—01 + 02 + 03, where the phases are defined via CkS* k = \ckSk\e l< f >k . Furthermore, 
we set k = |ci«i| |c2«2 1 IC3S3I . Note that the other offdiagonal elements follow from 
these and the fact that g is hermitian. 

From Eq. Q one can conclude two things. First, for a product states the 
absolute value of all the offdiagonal elements is the same. Second, the four phases 
are not arbitrary; they depend only on three parameters. 

To proceed, consider a linear functional like 

£(g,X) = $t(X l6ltS + X 2 g 2 ,7 + X 3 g 3fi + X 4 g 5A ), (8) 

where X = (Xi, X 2 , X 3 , X4) is a vector of complex coefficients and denotes 
the real part. The functional C is compatible with convex combinations of the 
quantum state, i.e., one has C[pg\ + (1 — p)g 2 ] = p£(Qi) + (1 — p)£(Q2)- For a 
pure separable state, C is given by 

£(|0)) = kF(X) with 

F(X) = K(Xi)cos(a + 6 + c) -%(X 1 )sm(a + b + c) + 3?(X 2 ) cos(a) 
-S(X 2 ) sin(a) + ^(X 3 ) cos(6) - 5S(X 3 ) sin(6) 
+$R(X 4 ) cos(c) - 3?(J*r 4 ) sin(c). (9) 

In order to obtain an extension of the separability condition in Eq. ^ we need 
two more facts. First, for given coefficients X one may compute the maximum 
given by 

C(X) =sup|J r (X)| (10) 

a,b,c 

For given values of X this can usually be computed analytically (see also below) 
or with a simple numerical optimization. 

Second, we have to characterize k for product states. For them we have 

« = \/£l,l£>4,4£>6,6£7,7 = \/£»2,2£>3,305,5£»8,8 

= ^Qk,kQ9-k,9-k for all 1 < k < 4. (11) 
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There are further equalities of this type (e.g., k = -^^2,2^3,3^4,4^5,5^6,6^7,7) and 
in view of Eq. ^ one might be tempted to use them. This, however, will not 
give stronger criteria: k will be used to deliver an upper bound, and in general 
a bound like x < \/af3j [corresponding to Eq. ((2J)] is weaker than the bound 
x < min{a, /3, 7} [corresponding to the last line in Eq. ([lip ]. In a similar way one 
can directly see that other possible bounds (as used in Ref. [13] ) can be derived 
from the terms in Eq. ([lip . 

Putting it all together, we can formulate: 

Observation. Let X be some coefficients and C{X) be defined as in Eq. H10\). 
Then, if g is fully separable the inequality 

\£{Q,X)\ < C(X) mm { ^/£>l,l£4,4£»6,6£>7,7, $Q2,2Q3,3Q5,5Q8,8, 

^Qk,kQ9-k,9-k {with 1 < k < 4)} (12) 

holds and violation of this inequality implies entanglement. 

To prove this criterion, note first that Eq. ()12p holds for pure product states. 
Furthermore, the left-hand side is a convex function of the state, while the right- 



hand side is concave 121 ] . Since a mixed fully separable state is a convex combi- 
nation of fully separable pure states, Eq. ([T2[) has to hold. 

It remains to discuss which parameters X should be chosen in order to de- 
tect a given state g. A simple choice is X ~ (g* 8 , g\ 7 , £3 6 , £5 4 ), since then the 
Cauchy-Schwarz inequality guarantees that C(g, X) is maximal for all X with 
the same normalization. The optimal choice, however, is to choose X such that 
\C(g, X)\/C(X) is maximal. 



4. Examples 

4-1. The bound entangled state from Eq. f3|) 

As a first example, let us discuss the state £>ak from Eq. ([3]). First, let us 
take a look at 

C[X = (5, a, /3, 7)] = sup[Jcos(a + b + c) + a cos(a) + (5 cos(6) + 7 cos(c)], (13) 

a,b,c 

where the coefficients 6, a, /3,7 are real. This case is important for GHZ-diagonal 
states, since for them all anti-diagonal elements are real. 

If these coefficients are all positive, it is clear that the maximum is C = \5\ + 
\ a \ + \P\ + I7I an d is attained at a = b = c = 0. The same value for C is obtained, 
if two or four of the coefficients are negative, because with transformations like 
a 1— > a + 7r one can flip two arbitrary signs of the cosine terms. 
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If one of the coefficients is negative (we choose 7 < for definiteness) one 
finds after some algebra that the optimum can be attained at 



sin(a) = =, cos(a) = ±yl — sin( 
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sin(6) = — — —=, cos(b) = — sin(6) 2 , 



sin(c) = cos(c) = ±-\A — sin(c) 2 , (14) 



with 



Q = — (a/35 + a/37 ~~ ~ P$l) (a/35 — a(3j + ac>7 — P5j) 
(a/35 — a/3~f — a5j + /35^)(a(35 + a/37 + a5j + /W7), 
i? = aPiSiaP - 5i)(aj - P8)(a6 - Pi). (15) 

This is the optimal solution, if the sin(...) are real and their absolute value is 
not larger than one. Furthermore, one has to distribute maximally one negative 
sign in the cos(...) terms such that the condition 5sin(a + b + c) + asin(a) = 
(originating from setting the derivative to zero) holds. If the conditions on sin(...) 
are not met, the optimum is just given by a choice of {a, b, c} G {0, ir} resulting 
inC = a + /3 + 7 + 5orC = a + /3 — 7 — 5 etc. The solution is of the same 
structure, if three of the coefficients a, P, 7, 5 are negative. 

Applying this to the state from Eq. ([3]) one finds that C[X = (1, 1, —1, 1)] = 
2\/2. According to Eq. Ifl2"j) . the state is entangled, if 8 > 2\/2a O a < 2y/2, 
which proves the numerical result. Recall that for larger values of d the state 



can be proven to be separable (see Ref. 16f] and also below), so the criterion is 
optimal. 

4-2. Random GHZ diagonal states 

As a second example, we consider randomly chosen GHZ diagonal states [see 
also Eqs. ([U [5])]. Here, we also want to estimate how good the criteria are, in 
the sense that we want to test whether some entangled states escape from the 
detection!! 

Consequently, we need a method to prove that a given state is separable. This 
is, in general, not straightforward, but there are two possibilities that we used: 
First, in Ref. [191 ] a simple iterative algorithm was proposed that can be used to 
prove different forms of separability for general quantum states. 

Second, for the special case of GHZ diagonal states one can also use the ideas 



proposed recently in Refs. [16|,ll7|. Let us explain them step by step. First, a state 



2 For results on certain subclasses of GHZ diagonal states see Refs. [lrl UtI [2^| . 
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like q ~ al + aZZ\ is separable, as it is diagonal in a product basis and has no 
negative eigenvalues. Then, consider an operator A = X2Z Zt + \3ZtZ '+ X4IZ Z '. 
The minimal eigenvalue of it is given by A_ = min{A2 + A3 + A4, A2 — A3 — 
A4, — A2 + A3 — A4, — A2 — A3 + A4} and A_ is never positive. Therefore, a state 
like q ~ |A_|1 + A is also separable, as it is diagonal in the computational basis 
and positive semidefinite. 

Considering the state in Eq. ([5]), this allows already to conclude that a GHZ 
diagonal state is separable, if |A_| + |As| + | Ae I + | A7 1 + 1 Ag | < 1. However, inspired 



by Ref. [17( we can go further: taking the operator B = X^XXX + XqYYX + 



XjYXY + X$XYY we may try to rewrite it as 

B=^[(A 1 +B 1 )®(A 2 +B 2 )®(A 3 +B 3 )+(A 1 -B 1 )®(A 2 -B 2 )®(A 3 -B 3 )], (16) 

with Ai = jUCOs(??i)X and B{ = fj,sin($i)Y. Then, the minimal eigenvalue of each 
of the terms in the right-hand side of Eq. (|16|) is given by — /i 3 , and the matrix 
g ~ /i 3 l + B is separable. 

If all Xi in B are positive, a solution of Eq. (|16p can be found with 



3/, , x \ \ \J (A5A6 + A 7 As)(AsA7 + A6As)(A5As + A6A7) 

V (A 5 , A 6 , A 7 , Ag) = ===== (17J 

V A5A6A7A8 

and the same solution holds if an even number of A« are negative. Therefore, in 
this case the state is separable, if |A_| + /i 3 < 1. For the case that an odd number 
of the Xi in B is negative, a decomposition as in Eq. (|16|) will lead to non- 
hermitian Ai and Bi which does not help. Then, only the separability condition 
|A_| + I A5 1 + I Ag I + I Av I + I Ag I < 1 remains. But the case that JliLs Aj < is 
exactly the case, for which it has been shown already in Ref. [16] that the PPT 
criterion is necessary and sufficient for full separability (note the different sign 
conventions in Ref. |l6l|). 

In some cases, these conditions can still be improved. For that, one may 
consider operators of the form 

X = p(Ai + B 1 )®(A2 + B 2 )®(A 3 + B 3 ) + (l-p)(A 1 -B 1 )®(A2-B 2 )® (i 3 - B 3 ) 

(18) 

with Ai = ctiX, Bi = fyY, Ai = otiX, Bi = fyY, and < p < 1. The minimal 

eigenvalues of the two terms are i] = — Yli(\J + 0f ) and fj = — Y\i(\J + /3f), 
respectively. Therefore, the state 

a K *~(p|77| + (i-p)|7)|)i + * (19) 

is separable. 

This state g sep is not diagonal in the GHZ basis as it contains also terms 
like XXY etc. But, by applying local operations one can make it GHZ diagonal 



Entangled 
[via Eq. CE2])] 


91.32 % 


NPT for some partition 


90.61 % 


PPT, but violating Eq. (Q2D 


0.71 % 


Separable 


8.68 % 


Via Eqs. (fl~6l[TT|) 


8.41 % 


Using in addition Eqs. (JTHJ US) 


0.27 % 


Total 


100 % 




100 % 



Table 1: Fractions of the randomly generated states, which are detected by the different criteria. 
See the text for details. 

without changing the weights of XXX, XYY, YXY, and YYX. In fact, after 
making the terms in Eq. [18] GHZ diagonal, a convex combination of two times 
Eq. [16] arises. 

Consequently, if we consider a given GHZ diagonal state g g h z = xl + B and 
find a separable state like g aep with the same weights [that is, paxotioij, + (1 — 
p)aia2&3 = A5 etc.] then g^ z must be separable if x > H 7 ?! + (1 — p)! 5 ?!]- The 
search for the appropriate X can easily be done numerically. This criterion is 
stronger than the one of Eqs. (fTUl [T7|) since it contains the latter as the special 
case p = 1/2. 

To test our criteria, we have generated 10 6 GHZ diagonal states, by choosing 
the eigenvalues pk in Eq. (j3]) randomly from the seven-dimensional simplex in 
ft 8 . Then, we tested the criterion of the positivity of the partial transposition 
[l|, 0] as well as our new criterion from Eq. (|12j) . For this criterion we found 
the optimal X via a simple numerical optimization. Note that for GHZ diagonal 
states the criterion in Eq. (|12p is strictly stronger than the PPT criterion: For 
instance, a state that is PPT with respect to the A\BC partition has to fulfill 
I f?i8 1 < y / i>44£55 which is a special case of Eq. (fT2"j) . For the states that were not 
detected by these criteria, we have tried to prove that they are separable using 
the ideas from above. The results are given in Table I. 

First, it is important that any state which was not detected by the new 
criterion in Eq. (|12|) was proven to be separable. This gives clear evidence that 
the criterion is a necessary and sufficient entanglement criterion for GHZ diagonal 
states. Moreover, the volume of the bound entangled states among the GHZ 
diagonal states can be estimated to be 0.7%. 

Finally, we also compared the general separability algorithm from Ref. (l9| 
with the special criteria for GHZ diagonal states in Eqs. (fT6HT9]l . It turns out that 
the general algorithm succeeds for 96% of the separable states and is therefore a 
powerful tool to prove separability, which can also be used if a state is not GHZ 
diagonal. 
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5. Bound entanglement in the vicinity of the W state 



The previous criterion was well suited for GHZ diagonal states. For other 
states, one may first apply local transformations in order to bring them close to 
a GHZ diagonal state. Nevertheless, this will not always succeed and one can 
therefore ask whether similar ideas can be used for other bound entangled states, 
which are not GHZ diagonal. To see that this is possible, let us consider the state 
investigated by P. Hyllus [241 ] . 
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(20) 



This state is close to the three-qubit W state \W) = (|001) + |010) + |100))/\/3, 
and is separable with respect to any bipartition, but not fully separable. The 
entanglement in this state is not detected by the methods presented in Ref. . 

To investigate the entanglement properties of these states, consider first filter 
operations of the form 



tot 



T ® T ® T with T 



i 

X 

x 2 



(21) 



Under this filtering operation the state transforms like ^ph(^) ^ ^totl^PH^j^tot 
~ Qpn(r,/x 6 ). Therefore, all states in this family share the same entanglement 
properties and one may focus on the case that r, = y/3/2, where the off-diagonal 
terms are maximal, 1 02,3 1 + 1 03,5 1 + 1 05,2 1 = 1/(1 + \/8/3) rj 0.38. 

In order to derive a separability criterion, we consider a pure product state as 
in Eq. (jSJ) with | si | 2 1 S2 1 2 1 S3 1 2 = £s,8 < £■ Then, for some i G {1,2,3} the bound 
I I < \fs must hold, and from that it follows after a short calculation (using 
\cj\ < 1) that 1 ,02, 3 1 + I £3,5 1 + 1 05,2 1 — \/£ + 1/4. Therefore, any pure separable 
state obeys 



I i»2,3 1 + I 03,5 1 + 1 05,2 1 



■+1/4. 



(22) 



The left-hand side of this inequality is convex, while the right-hand side is con- 
cave, so this inequality holds also for mixed separable states. The state £ph with 
r] = y/3/2 clearly violates it (also if some noise is added), so the states in the 
family gpu(f]) must be entangled. 
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6. Conclusion 

In conclusion, we have presented a criterion for full separability in terms of 
inequalities for the density matrix elements. The criterion seems to be necessary 
and sufficient for GHZ diagonal three-qubit states and allows to characterize 
bound entangled states. Finally, we showed that our ideas can also be extended 
to different families of bound entangled states close to the W state. 

There are several possibilities to extend our work in the future. First, an 
analytical proof of the necessity of the criterion for GHZ diagonal states would 
be desireable. Then, an explicit investigation of the criteria for four or more 
qubits would be of interest. Finally, one could try to extend our analysis to the 
more general class of states which are diagonal in a graph state basis. 

We thank B. Jungnitsch, P. Hyllus, A. Kay, M. Kleinmann, T. Moroder, S. 
Niekamp and M. Seevinck for discussions. This work has been supported by the 
FWF (START Prize and SFB FOQUS) and the EU (NAMEQUAM). 
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